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ABSTRACT 

We investigate the evolution of the galaxy stellar mass function at high-redshift {z > 5) 

using a pair of large cosmological hydrodynamical simulations: MassiveBlack and 
MassiveBlack-II. By combining these simulations we can stiidy the properties of galax- 
ies with stellar masses greater than 10^ Mq/;,^^ and (co-moving) number densities of 
logiQ(^ [Mpc~^ dex~^ /i^]) > —8. Observational determinations of the galaxy stellar 
mass function at very-high redshift typically assume a relation between the observed 
UV luminosity and stellar mass-to-light ratio which is applied to high-redshift samples 
in order to estimate stellar masses. This relation can also be measured from the sim- 
ulations. We do this, finding two significant differences with the usual observational 
assumption: it evolves strongly with redshift and has a different shape. Using this 
relation to make a consistent comparison between galaxy stellar mass functions we 
find that at z = 6 and above the simulation predictions are in good agreement with 
observed data over the whole mass range. Without using the correct UV luminosity 
and stellar mass-to-light ratio, the discrepancy would be up to two orders of magni- 
tude for large galaxies > 10^° Mq h~^. At 2; = 5, however the stellar mass function for 
low mass < 10^ galaxies is overpredicted by factors of a few, consistent with 

the behaviour of the UV luminosity function, and perhaps a sign that feedback in the 
simulation is not efficient enough for these galaxies. 

Key words: galaxies: evolution - galaxies: formation - galaxies: starburst - galaxies: 
high-redshift - ultraviolet: galaxies 



1 INTRODUCTION 

The observational exploration of the high-redshift [z > 2) 
Universe has been driven, over the past 10-15 years, predom- 
inantly by deep Hubble Space Telescope (HST) surveys. Deep 
Advanced Camera for Surveys (ACS) observations alone (of 
the HUDF for example) permitted the identification of large 
numbers of galaxies at z = 2 — 6 (e.g. Bunker et al. 2004, 
Beckwith et al. 2006, Bouwens et al. 2007). While some 
galaxies aX z > 7 were identified using ACS and near-IR 
Camera and Multi-Object Spectrometer (NICMOS) obser- 
vations (e.g. Bouwens et al. 2008) or ground based imaging 
(e.g. Bouwens et al. 2008, Ouchi et al. 2009, Hickey et al. 
2010) the very-high redsiiift Universe was only truly opened 
up by the installation of Wide Field Camera 3 (WFC3) in 
2009. WFC3 near-IR (1.0 - 1.6/im) observations allow the 
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identification of star forming galaxies to z — 7—8 (e.g. 
Bouwens et al. 2010, Ocsch et al. 2010, Bunker et al. 2010, 
Wilkins et al. 2010, Wilkins et al. 2011a, Lorenzoni et al. 
2011, Bouwens et al. 2011b) and potentially even to « ~ 10 
(Bouwens et al. 2011a, Oesch et al. 2012). 

By combining ACS optical and NICMOS or WFC3 
near-IR imaging with Spitzer IRAC observations it becomes 
possible to probe the rest-frame UV-optical spectral energy 
distributions of galaxies at 2 = 4 — 8 (e.g. Eyles et al. 2005, 
Gonzalez et al. 2012). Rest-frame optical photometry is cru- 
cial to accurately determine stellar masses (e.g. Eylcs et 
al. 2007, Stark et al. 2009, Labbe et al. 2010, Gonzalez et 
al. 2011). With a sufficiently large, well defined sample of 
galaxies it is possible to study the galaxy stellar mass demo- 
graphics, and in particular the galaxy stellar mass function 
(e.g. Gonzalez et al. 2011). The galaxy stellar mass function 
(GSMF) is a fundamental description of the galaxy popu- 
lation and is defined as the number density of galaxies per 
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Table 1. Main characterics of Massive Black and Massive Black- 
II simulations. Both simulations included dark matter, SPH, a 
multiphase model for star formation, and a model for black hole 
accretion and feedback. The number of particles Npa^t is given, 
the size of the simulation box L^oxi the gravitational softening 
length e, the number of cores used Ncores and the final redshift 
Zf. Both runs were started at 2: = 159 and used 6 threads/MPI 
task. For MassiveBlack-II the number of cores and threads used 
was optimized as it progressed. 



Run Npart Lbox e 

(Mpc/h) (kpc/h) 



MassiveBlack 2 X 3200^ 533 5.0 4.75 

MassiveBlack-II 2 X 1792=* 100 1.85 



major improvement over previous versions of GADGET is in 
the use of threads in both the gravity and SPH part of the 
code which allows the effective use of multi core processors 
combined with an optimum number of MPI task per node. 
The MassiveBlack simulation contains Npart = 2 x 3200^ — 
65.5 billion particles in a volume of 533Mpc//i on a side 
with a gravitational smoothing length e = 5.0kpc//i in co- 
moving units. The gas and dark matter particle masses are 
TTig = 5.7 X 10' Mq and ttidm = 2.8 x 10* respectively. 
The simulation has currently been run from 2 = 159 to 
z = 4.75 (beyond our original target redshift of z — 6). For 
this massive calculation it is currently prohibitive to push 
it to 2 = as this would require an unreasonable amount 
of computational time on the world's current fastest super- 
computers. The simulated redshift range probes early struc- 
ture formation and the emergence of the first galaxies and 
quasars. 

MassiveBlack- 11 (see Khandai et al. in-prep for an 
overview) is a smaller volume but the mass and spatial res- 
olution are better than MassiveBlack by a factor of 25 and 
2.7 respectively. The smaller volume means that a smaller 
part of the high mass function can be sampled and that in 
the mass range where it overlaps with MassiveBlack it can 
be used to check for convergence as well as to extend our 
predictions towards the low mass end. This is the largest 
volume ever run at this resolution with a final redshift of 
z = 0. 

These runs contain gravity and hydrodynamics but 
also extra physics (subgrid modeling) for star formation 
(Springel & Hernquist 2003), black holes and associated 
feedback processes (Di Matteo et al. 2008, Di Matteo et 
al. 2012). The cosmological parameters used were; the am- 
plitude of mass fiuctuations, ag, — 0.8, spectral index, 
ris = 0.96, cosmological constant parameter Qa ~ 0.74, 
mass density parameter Urn = 0.26 , baryon density pa- 
rameter Qb = 0.044 and h = 0.72 (Ifubble's constant in 
units of 100km s"^Mpc"^ WMAP5) for MassiveBlack. For 
MassiveBlack-II we instead used f^A ~ 0.725, and fim ~ 
0.275 (according to WMPA7). 

Catalogues of galaxies are made from the simulation 
outputs by first using a friends-of-friends groupfinder and 
then applying the SUBFIND algorithm (Springel 2001) to 
find gravitationally bound subhalos. The stellar component 
of each subhalo consists of a number of star particles, each 
labelled with a mass and the redshift at which the star par- 
ticle was created. 

To generate the spectral energy distribution (SED), 
and thus broad-band photometry, of each galaxy we sum 
the SEDs of each star particle (weighted by the particle 
mass). The SED of each star particle is generated using the 
Pegase.2 stellar population synthesis (SPS) code (Fioc & 
Rocca-Volmerange 1997,1999) taking account of their ages 
and metallicities. Nebula (contiimum and line) emission is 
also added to each star particle SED, though this has a neg- 
ligible effect on the UV photometry considered in this work. 
In addition we apply a correction for absorption in the in- 
tergalactic medium (IGM) using the standard Madau et al. 
(1995) prescription (though again this has a negligible effect 
on this work) . Throughout this work we measure the broad- 
band UV luminosity using an idealised rest-frame top-hat 
filter at A = 1500 ± 200A. A rest-frame filter is chosen to 
allow a consistent comparison between samples at different 



logarithmic stellar mass bin. The first moment of the GSMF 
corresponds to the cosmic stellar mass density. 

Here we use state-of-the-art cosmological hydrodynam- 
ical simulations of structure formation {MassiveBlack and 
MassiveBlack-II) to investigate their predictions of the 
GSMF and compare it with current constraints. These runs 
are large, high resolution simulations, with more than 65.5 
billion resolution elements used in a box of roughly cu- 
bic gigaparsec scales (for MassiveBlack), making it by far 
the largest cosmological Smooth Particle Hydrodynamics 
(SPH) simulation to date with full physics of galaxy for- 
mation (meaning here an inclusion of radiative cooling, 
star formation, black hole growth and associated feedback 
physics) ever carried out. The combination of the two sim- 
ulations allows us to probe galaxies with stellar masses 
greater than 10* M0 and (co-moving) number densities 
of logio(0[Mpc-^dex-i/i3]) > -8, a range well matched 
with current observations at high-redshift. 

This article is organised as follows: in Section [2] we in- 
troduce the MassiveBlack and MassiveBlack-II simulations. 
In Section[3]we explore the predicted evolution of the galaxy 
stellar mass function, how both the intrinsic and observed 
luminosities correlate with the stellar mass-to-light ratio and 
in §3.5| compare galaxy stellar mass functions to recent ob- 
servations. Finally, in Section |4] we present our conclusions. 

Throughout this work magnitudes are calculated using 
the AB system (Oke & Gunn 1983). We assume Salpeter 
(1955) stellar initial mass function (IMF), i.e.: ^(m) — 
dN/dm oc m-2■^^ 

2 MASSIVEBLACK AND MASSIVEBLACK-II 

2.1 Simulation runs: Massive Black and Massive 
Black-II 

Our new simulations (see Table [l] for the parameters of 
the simulation) have been performed with the cosmological 
TreePM-Smooth Particle Hydrodynamics code P-GADGET, 
a hybrid version of the parallel code GADGET2 (Springel 
2005) which has been extensively modified and upgraded to 
run on the new generation of Petaflop scale supercomputers 
(e.g. machines like the upcoming Blue Waters at NCSA). The 
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redshifts. The shape of this filter is selected for convenience, 
but closely reflects the profile of near-IR bandpasses which 
are available to measure the rest-frame UV flux at high- 
redshift. 

We note that our work is complementary to the recent 
simulation predictions of the galaxy stellar mass functions 
of Jaacks et al. (2012), who compare results for a suite of 
smaller simulations to the Gonzalez et al. (2011, hereafter 
Gil) observational data. Our work differs in extending to a 
lower redshift, correcting for the effect of an evolving ratio 
of UV luminosity to mass to light ratio, and also for the 
inclusion of supermassive black hole formation and feedback 
in our simulations. We discuss the Jaacks et al. (2012) results 
further below. 



3 THE GALAXY STELLAR MASS FUNCTION 

Measuring the GSMF from outputs of the MassiveBlack and 
MassiveBlack-II simulations is straightforward, given that 
the total masses of star particles in each galaxy are known. 
Before making a comparison to observational data, however, 
we must remember that observed UV luminosities were used 
(e.g. by Gonzalez et al. 2011, hereafter Gil) to compute the 
published observed GSMFs. This means examining the re- 
lationship between UV luminosity and stellar mass to light 
ratio in the simulation and using this information in our 
comparison to observations. In this section, we do this, af- 
ter first presenting the GSMF measured directly from the 
simulations. 
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Figure 1. The galaxy stellar mass function measured from the 
MassiveBlack (dashed lines) and MassiveBlack- II (solid lines) 
simulations for z £ {5, 6, 7, 8, 9, 10}. The two right hand axes show 
the number of galaxies in the MassiveBlack and MassiveBlack-II 
volumes. 



3.1 Galaxy Stellar Mass Function from 
simulations 

The evolution of the > 10* Mq galaxy stellar mass func- 
tion from 2 = 10 — >■ 5 predicted by MassiveBlack and 
MassiveBlack-II is shown in Fig. [l] The shape of the sim- 
ulated GSMF is a declining distribution with mass and, at 
least at 2; = 5, exhibits a sharp cut off at high-masses. Val- 
ues of the number density (j> are also tabulated in Table [2] in 
various logarithmic mass intervals. 

Figure [T] also demonstrates the evolution in the normal- 
isation of the GSMF. At 2 = 10 there are only ~ 500 galax- 
ies with stellar masses > 10* Mq in the MassiveBlack- 
II volume (10*' Mpc^ /i"^), while at 2 = 5 this has in- 
creased to ~ 135, 000 (x270). The shape of the GSMF also 
evolves strongly; while the number of galaxies with masses 
> 10* M© h'^ increases by a factor of x 270 from 2 = 10 — i> 5 
the number of galaxies with masses > 10^" Mq h^^ increases 
by a factor of x5000. 

The evolution of the simulated galaxy stellar mass func- 
tion is stronger than that exhibited by the UV luminosity 
function. This reflects the fact the average UV mass-to-light 
ratio of galaxies also increases 2 = 10 — >■ 5 (as demonstrated 
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3.2 



Observational Estimation of the Galaxy 
Stellar Mass Function 



spectral energy distributions of high-redshift galaxies. Rest- 
frame optical photometry is vital to determine accurate stel- 
lar masses. Several studies have recently attempted to mea- 
sure the stellar masses of high-redshift Lyman-break selected 
galaxies (e.g. Eyles et al. 2007, Stark et al. 2009, Labbe et al. 
2010, Gonzalez et al. 2011). With a sufhciently large sample 
and a handle on the incompleteness issues it is also possible 
to study the galaxy stellar mass function (e.g. Stark et al. 
2009, Labbe et al. 2010, Gonzalez et al. 2011). 

To understand how to make simulation predictions 
it is useful to examine exactly how the Gil GSMF is 
constructed. The Gil study draws a sample of galaxies 
from the observed UV luminosity functions (LFs) at 2; € 
{3.8,5.0,5.9,6.8} (using Bouwens et al. 2007, 2011). These 
UV luminosities are converted into stellar masses using 
the observed UV Luminosity (I/1500, 06s) - stellar mass-to- 
light ratio (Af/Li5oo,o()s) distribution measured at 2 ~ 4. 
This relation is fairly well fit by a power lawQ such that 
Af/I/1500, 06s oc L°'^ (i.e. the stellar mass-to-light ratio in- 
creases with observed UV luminosity) . While this relation is 
calibrated at 2 = 4 Gil note that that it appears to fit ob- 
servations of stellar masses and luminosities at 2 G {5,5.9}. 
However at these redshifts the sample sizes are small (78 
and 28 galaxies at 2 ~ 5 and 2; ~ 6 respectively) and there 
is a large degree of scatter. 



By combining HST optical and near-IR observations (from 
ACS and NICMOS or WFC3) with Spitzer IRAC photom- 
etry it is possible to measure the rest-frame UV-optical 



^ Though the power law fit is not used to determine the GSMF. 
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Table 2. The number density (in units of Mpc ^ dex ^ h^) of galaxies in various logarithmic mass intervals ([9.5,10.0) = 9.5 < 
logi(,(J\/) < 10.0, where M has units Mq h'^) for z e {5,6,7,8,9, 10}. Where there are no objects within the mass interval the number 
density is replaced by an upper limit corresponding to n < 1 (i.e. <)> < 1/V). 



Mass Interval 




lo 


gio{</'[Mpc- 


-3dex-i h^]) 






logio([M0ft-i]) 


z = 5 


z = 6 


z = 7 


z = 8 


2 = 9 


2 


= 10 


MassiveBlack Volume = (533 Mpc /i-i)^ 


[9.5, 10.0) 


-2.85 


-3.50 


-4.25 


-5.13 


-6.07 




-7.88 


[10.0,10.5) 


-3.69 


-4.46 


-5.35 


-6.43 


-7.88 


< 


-8.18 


[10.5,11.0) 


-4.55 


-5.46 


-6.80 


-7.88 


< -8.18 


< 


-8.18 


[11.0, 11.5) 


-6.23 


-7.88 


< -8.18 


< -8.18 


< -8.18 


< 


-8.18 


[11.5, 12.0) 


< -8.18 


< -8.18 


< -8.18 


< -8.18 


< -8.18 


< 


-8.18 


MassiveBlack-II Volume = (100 Mpc /i~i)3 


[8.0,8.5) 


-0.70 


-1.06 


-1.46 


-1.92 


-2.44 




-3.03 


[8.5,9.0) 


-1.27 


-1.69 


-2.15 


-2.69 


-3.31 




-4.06 


[9.0,9.5) 


-1.95 


-2.42 


-3.01 


-3.71 


-4.47 




-5.10 


[9.5, 10.0) 


-2.76 


-3.37 


-4.12 


-4.80 


-5.70 


< 


-6.00 


[10.0,10.5) 


-3.57 


-4.27 


-4.74 


-5.70 


< -6.00 


< 


-6.00 


[10.5, 11.0) 


-4.40 


< -6.00 


< -6.00 


< -6.00 


< -6.00 


< 


-6.00 



3.3 The relation between UV luminosity and the 
stellar mass-to-light ratio in simulations 

As noted above, the Gil study uses the distribution of stel- 
lar masses and UV luminosities measured at z ~ 4 to effec- 
tively convert the observed UV luminosity function into a 
galaxy stellar mass function. To make a proper simulation 
prediction we must take into account any difference between 
the relation between UV luminosity and the stellar mass-to- 
light ratio used by Gil and that in the simulations. 

Figure [2] shows the relationship between the intrinsic 
UV luminosity (Lisoo) and mass-to-light ratio (A//L1500) 
at z e {5, 6, 7, 8, 9, 10} predicted by MasstveBlack-II. This 
relationship is (over the full mass range) approximately flat 
(i.e. the intrinsic stellar mass-to-light ratio is constant) and 
is significantly different from the i\f/-£/i5oo, 06s oc L"'^ relation 
found by Gil. Jaacks et al. (2012) plotted the rest frame 
UV magnitude against stellar mass in their simulations, also 
finding a fiatter relationship that than used by Gil. The 
lower-panel of Fig. [2] shows that the relationship between 
the intrinsic UV luminosity and stellar mass-to-light ratio 
also varies strongly with redshift, increasing by 0.6 dex from 
z = 10 5. 

It is also interesting to note from Figure [2] that it ap- 
pears the intrinsic UV luminosity of galaxies with Lisoo > 
IQ^^ ergs~^ can alone be used to estimate the stellar 
mass with an accuracy of ~ 50%. This contrasts sharply 
with the low-redshift Universe where star formation has ter- 
minated in many systems (particularly massive ellipticals) 
rendering the UV luminosity to be negligible. The strong 
correlation between UV luminosity and stellar mass reflects 
the fact that virtually all galaxies at high-redshift (in the 
MassiveBlack and MasstveBlack-II simulations) continue to 
actively form stars. 



nosities (as used in Fig. [2|. Attenuation due to dust both 
decreases the UV luminosity (i.e. I/1500, 06s < iisoo) and 
increases the stellar mass-to-light ratio (i.e. Af/Lisoo, ots > 
M/Lisoo) relative to their intrinsic values. A positive corre- 
lation between luminosity and dust attenuation would then 
introduce a positive correlation between M /Li5oo,obs and 
the observed UV luminosity. 

The measurement of dust attenuation at high-redshift is 
challenging. Far-IR observations, and optical emission lines, 
are generally inaccesible for the bulk of the galaxy popula- 
tion at high-redshift leaving only the UV continuum slope /3 
as a diagnostic (e.g. Meurer et al. 1999, Wilkins et al. 2012a, 
Wilkins et al. submitted). A number of recent studies have 
attempted to constrain the relationship between j3 and the 
observed UV luminosity at high-redshift though with some 
conflicting results (e.g. Stanway, McMahon, & Bunker 2005, 
Bouwens et al. 2009, Wilkins et al. 2011b, Dunlop et al. 2012, 
Bouwens et al. 2012, Finkelstein et al. 2012). Bouwens et al. 
(2009), Wilkins et al. (2011b) and Bouwens et al. (2012) an 
increase in /? with observed luminosity. Dunlop et al. (2012) 
and Finkelstein et al. (2012) on the other hand found little 
evidence of variation of /3 with luminosity (see Wilkins et al. 
submitted for a detailed comparison). 

Adopting the relationship(s|^between /3 and luminosity 
found by Bouwens et al. (2012) and utilising the Meurer et 
al. (1999) calibration (between the observed UV continuum 
slope /3 and UV attenuation) we can determine the relation- 
ships between the observed UV luminosity (iisoo.obs) and 
observed mass-to-light ratio at z € {5, 6, 7} as predicted by 
MassiveBlack and MassiveBlack-II. These are shown in Fig. 
[3] The most significant change (relative to that found for the 
intrinsic luminosities and mass-to- light ratios) is that the re- 
lationship between -L1500, 06s and Af/Li5oo,obs is no longer ap- 



3.4 The effect of dust attenuation 

The Gil relation is however based on the observed (i.e. dust 
attenuated luminosities) as opposed to the intrinsic lumi- 



^ If a luminosity invariant dust correction was assumed the shape 
of the observed U V luminosity - mass-to-light ratio relation would 
remain the same (though the average observed mass-to-light ratio 
would increase). 
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Figure 2. The relationship between the intrinsic UV luminosity 
and stellar mass-to-light ratio at z S {5,6,7,8,9,10} predicted 
from MassiveBlack-II. In both panels the points denote the me- 
dian value of the mass-to-light ratio in each luminosity bin. In the 
upper-panel the 2-d histogram shows the density of sources on a 
linear scale and the error bars show the range encompassing the 
central 68.2% of galaxies. The arrow in the upper-panel shows the 
effect of dust attenuation (the labels denote values of A1500). The 
dashed line in both panels shows M/L1500 on L" '^ which provides 
a good fit to the distribution used by Gonzalez et al. (2011) to 
determine stellar masses from observed UV luminosities. 

proximately constant but is instead strongly positively cor- 
related, at least at Mi5oo,obs < —19.5 . At Afisoo.obs < —19.5 
the slope of this relation is 7 = 0.5 — 0.8 (where 7 is defined 
such that M/Z/i5oo,o6s oc L^) (c.f. 7 = 0.7 found by Gil at 
z = 4). This suggests the physical cause of the strong ob- 
served correlation between UV luminosity and mass-to-light 
ratio is caused almost solely by the correlation of dust at- 
tenuation with luminosity. At lower-luminosities the relation 
flattens (7 < 0.2). This arises due to the diminishing effect 
of dust at lower- luminosities, i.e. the I/i5oo,o6s - M / Lisoo.ota 
begins to reflect the (virtually fiat) intrinsic relation. 

3.5 Comparison with observations 

We are now in position to compare the MassiveBlack and 
MassiveBlack-II results to observations. We follow a proce- 
dure similar to Gil but using the simulated relation between 
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Figure 3. The relationship between the dust attenuated (ob- 
served) UV luminosity and stellar mass-to-light ratio at 2 £ 
{5, 6, 7} predicted from MassiveBlack-II aad using Bouwens et al. 
(2012) to relate dust attenuation to the observed UV luminosity. 
The points denote the median value of the mass-to-light ratio in 
each bin while the vertical error bars (at 2 = 5) denote the 68.2% 
confidence interval. The diagonal lines denote M/L1500 obs <^ L'' 
for 7 = {0.1, 0.2, 1.0}. The dashed line denotes 7 = 0.7. 



UV luminosity and mass-to-light ratio. We construct a vol- 
ume limited sample (referred to below as "B07/B11-I-MB 
MTOL") of galaxy UV luminosities using the Bouwens 
et al. (2007, 2011) observed UV luminosity functions. We 
then convert the observed UV luminosity of each galaxy to 
a stellar mass using the relation between luminosity and 
stellar mass-to-light ratio (M/Li5oo,ot.s) predicted by the 
MassiveBlack-II simulation (combined with the empirical 
dust correction described above) and construct a galaxy 
stellar mass function. These galaxy stellar mass functions 
are shown at z G (5, 6, 7} in Fig. [4] We also show in Fig. [4] 
the GSMFs predicted by MassiveBlack/ MassiveBlack-II axiA 
those determined by Gonzalez et al. (2011, hereafter Gil) 
at z G (5, 6, 7} (the z ~ 7 GSMF comes from Labbe et al. 
2010 but is also presented in Gil). 

From an examination of Fig. |4] it is clear that the 
B07/Bll-fMB MTOL sample shows a much closer corre- 
spondence to the simulations compared to Gil. This es- 
sentially reflects the good overall agreement between the 
simulated UV LF and the observations, at least at high- 
luminosities. The flattening of the relation between L1500, 06s 
and the mass-to-light ratio at low-luminosities does go some 
way to explaining the difference in the shape of the sim- 
ulated and observed galaxy stellar mass functions. More 
importantly however is the strong redshift evolution: from 
z = 5 — >■ 7 the calibration relating the the observed UV lu- 
minosity to the mass-to-light ratio decreases by 0.3 — 0.5dex 
(depending on the luminosity). Because the Gil study as- 
sumed no redshift evolution (instead utilising a calibration 
based on observations at z ~ 4 to convert UV luminosities 
to stellar masses at z = 4 — 7) this would cause the stellar 
masses to be overestimated. Because the GSMF declines to 
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ever recent simulations of isolated galaxies (e.g. Hopkins, 
Quataert & Murray 2012) have shown that, in the presence 
of feedback, restricting star formation to molecular gas or 
modifying the cooling function has very little effects on the 
star formation rates. By contrast changing feedback mecha- 
nism or associated efficiencies translates in large differences 
in final stellar mass densities. Based on these recent results 
(albeit on idealized simulations) we are prone to interpret 
our discrepancy at the low mass end to details in the stel- 
lar feedback model (and in particular to its efficiency which 
may be too low). 

Comparing to the simulation results of Jaacks et al. 
(2012) (which do not include AGN modelling), we see that 
a similar sign to the disagreement with observations at low 
mass. At the high mass end, we have shown that correcting 
for the UV luminosity-mass to light ratio assumed brings 
the observations and simulations into agreement, and this 
would also be likely to work for the Jaacks et al. results. 
Finally, it is also worth noting that the Gil GSMF evolves 
only very mildly from 2: = 5 — )■ 7. Indeed, the stellar mass 
density (which is the first moment of the GSMF) of galaxies 
with > 10* M0 is virtually fiat z = 5 — ^ 7. This is surpris- 
ing given that all the galaxies contributing to the GSMF at 
these redshifts/masses are likely actively forming stars (by 
virtue of being UV selected) and suggests either the high- 
redshift GSMF is overestimated or the lower-redshift GSMF 
underestimated. 



Figure 4. The galaxy stellar mass function predicted by Massive- 
Black (dashed lines) and MassiveBlack-II (solid lines) compared 
with observations at 2 £ {5, 6, 7} (top, middle, and bottom panels 
respectively). The open symbols in each panel show the predic- 
tion for the GSMF using the Bouwens et al. (2007, 2011) observed 
UV LP and a relationship between stellar mass and luminosity 
derived from MassiveBlack-II. The filled grey points show the 
GSMF from Gonzalez et al. (2011), which was estimated using 
a non-evolving relationship between UV luminosities and stellar 
mass to light ratios. Note that the units now implicitly assume 
h = 0.7. 



high-masses this would cause the number density of sources 
at any mass to overestimated. 

We also note from Fig. |4] that at 2: = 5 (and to a 
lesser extent at z = 6) this process does not fully reconcile 
the GSMF at low-masses. At 2: = 5 MassiveBlack-II over- 
predicts the faint-end of the UV luminosity function rela- 
tive to the observations of Bouwens et al. (2007) by around 
a factor x5 at Afisoo ~ —18. This is difficult to reconcile 
observationally without requiring the application of a much 
larger completeness correction. It therefore suggests that the 
discrepancy has its roots in the MB/ MB-II modelling as- 
sumptions. This disagreement occurs in low mass galaxies 
which are much less affected by AGN feedback and hence 
more sensitive to the details of the star formation model 
and stellar feedback. For example our model does not in- 
clude any treatment of the molecular gas component such 
as in e.g. Krumholtz and Gnedin (2011) which would tend to 
suppress star formation rates in lower mass galaxies. How- 



4 CONCLUSIONS 

We have investigated the high-redshift {z = 5—10) evolution 
of the galaxy stellar mass function (GSMF) using a pair of 
large cosmological hydrodynamic simulations MassiveBlack 
and MassiveBlack-II. Over the redshift range 2 = 10 — >■ 5 we 
find both the normalisation and shape of the GSMF evolves 
strongly with the number density of massive galaxies (> 
10® Mq) increasing by a factor of around x300. 

By combining Hubble Space Telescope optical and near- 
IR observations (from ACS, NICMOS and WFC3) with 
near-IR IRAC photometry from the Spitzer Space Telescope 
it is possible to identify and measure the stellar masses of 
galaxies at very-high redshift, and thus constrain the GSMF 
(e.g. Gonzalez et al. 2011). While the simulated GSMF at 
z — 5 provides reasonable agreement with the Gonzalez et 
al. (2011) observations at > 10®'^ Mq, at low-masses and at 
2 > 5 there is a significant discrepancy. The disagrement at 
low-masses at 2 = 5 is also refiected in the UV luminosity 
function (LF) at low-luminosities. However, at 2 > 5 the dis- 
crepancy appears to arise due to a difference in the assumed 
relationship between the observed UV luminosity and mass- 
to-light ratio. Gonzalez et al. (2011) applies a relationship 
calibrated at 2 ~ 4, however we find that the relation, while 
having a similar form (i.e. that the mass-to-light ratio is pos- 
itively correlated with the observed UV luminosity), evolves 
strongly with redshift. Applying a calibration based on the 
simulated distribution of UV luminosities and stellar masses 
to the observed UV luminosity functions yields galaxy stel- 
lar mass functions which closely reflect those predicted by 
the simulations. This simply reflects the good agreement be- 
tween the observed and simulated intrinsic UV luminosity 
functions. 
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